Vacuum Energy and Casimir Force in a Presence of Skin-depth 

Dependent Boundary Condition 

S.L. Lebedev Q 

The vacuum energy-momentum tensor (EMT) and the vacuum energy corresponding to 
massive scalar field on 3^ x [0, 1] x 'Sl D ~ 2 space-time with boundary condition involving 
a dimensional parameter (<5) are found. The dependent on the cavity size I Casimir en- 
ergy Ec is a uniquely determinable function of mass m, size I and "skin-depth" 5. This 
energy includes the "bulk" and the surface (potential energy) contributions. The latter 
dominates when I ~ 5. Taking the surface potential energy into account is crucial for the 
coincidence between the derivative —dEc/dl and the //-component of the vacuum EMT. 
CN ' Casimir energy Ec and the bulk contribution to it are interconnected through Legendre 

5-h ■ transformation, in which the quantity <5 _1 is conjugate to the vacuum surface energy mul- 

tiplied by 5. The surface singularities of the vacuum EMT do not depend on I and, for 
even D, 5 = or oo, possess finite interpretation. The corresponding vacuum energy is 
' finite and retains known analytical dependence on the dimension D. 

(N ; 

j> . 1. Introduction 

a\ ; 

The determination of vacuum energy of fields confined to finite volumes is the basic part in 
^ the calculations of Casimir force as well as of some bag characteristics of hadrons |], • At the 

same time the renormalization approach appears to be not unique: the counterterms necessary 
to make the vacuum energy finite can depend on the cavity size so that the regularization 
dependence of Casimir energy appears || [§]. Concrete calculations for the scalar case and 
Dirichlet or Neumann boundary condition on a spherical boundary have been performed in the 
Qh! recent papers || [7|. General consideration || shows however that at least for these bound- 

ary conditions the surface singularities of the vacuum EMT (and, hence, the aforementioned 
counterterms) can not depend on the size of the cavity. The detection of such dependences 
in the papers listed above is a consequence of a coincidence between the local (curvature) and 
global (size) parameters determining the spherical boundary. 

We consider the vacuum characteristics of massive scalar field defined on the domain 5i t x 
[0, /] x ^t D ~ 2 . In the absence of curvature we introduce "skin-depth" parameter 5 in boundary 
condition 

di(p(t, 0, x ± ) = -<p(t, 0, x ± ), (p(t, I, x ± ) = (1) 
This boundary condition implies the presence of the surface potential energy |H| 

l|^- 2 x ± ^ 2 (t,0,x ± ) (2) 

in the total hamiltonian of the model. The latter is, as usual, space integral of 00-component 
of the energy-momentum tensor 

= - [d^, d v (p] + - ^g^idif) 2 + ^g^ u m 2 (p 2 + £d x {g^ [</?, <9 A y?] + - [if, d v tp} + ) (3) 
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2 )We use the system of units where h = c = 1. The signature of metric is (+, .). x M = (t,xi,x±). 

Brackets [,]+ mean anticommutator. 
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where the last divergence term should be taken with coefficient £ = 1/4. This is the only value 
of £ guaranteeing the diagonality of hamiltonian and its conservation JTOl . 



2. Energy-momentum tensor and energy of vacuum 

The normalized solutions of Klein-Gordon-Fock equation have the form: 

<Pk{x) = (2co>)~2 (2-k) 1 ~~ exp [—iut + iqx.±]ipk(xi) , 



where to = y/m 2 + q 2 + k 2 , and the descrete set of functions 



i) k {x) = N k sin k(x -I), N k 



I i 



sin 2kf 
2kl 



-1/2 



(4) 



(5) 



represents the resonator modes. Wave number k = z/8 (1 = 1/8) should be determined from 
the spectral equation 

A(z) = z _1 tgzr+l = 0, (6) 

following from boundary condition (1). The basis of functions (4) is used to define creation- 
annihilation operators and vacuum state in a standard manner |T|, [TTJ 

With the help of Cauchy theorem applied to meromorphic function A(z), it is possible 
to transform the sums over transcendential roots of eq.(6) into corresponding integrals. The 
scheme of calculation fits in the formula 



E 

z nS >0 



f(z n s) 



I 



sin 2z n sl 



f(z) dz 



I 7 (i - *)[/(«) - /(-«)] 



2vri 



shtl + t ch tl 



dt, 



(7) 



where f(z) is a function analitical in the right half-plane, and z n s (n = 1, 2, • ■ ■) are the (real) 
roots of eq.(6). With the use of eq.(7) one could attach the integral representation to the un- 
renormalized vacuum EMT. The renormalization is reduced to the subtraction of Minkowskian 
space contribution (D-regularization [lTSJ supposed): 



y?(4vr) Z) / 2 Jo 



POO 

/ dk M D ~ 3 diag 
Jo 



—M 2 
k 2 (l-D) 
M 2 



M=(m 2 + k 2 ) 1 ' 2 . The final answer i.e. renormalized EMT of vacuum, looks like 

(1 - t)e- tT 



5.1 
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dt(t 



shtl + t ch tl 



diag 



ch 2tx' - 1 
t 2 (l~D) 



t 2 -n 2 
- ch2tx' 



+ e 



-2tx' 



diag 



1 




(9). 

(D — 2)-dimensional diagonal of {T^ u )s t i corresponding to xi contains expressions of the same 
type. The following notations were used in eq.(9): 



K 



8 



l-D 



D 



2(4^(1^) 



\i = m8, x' = (I — xi)/8. 



(10) 



3 ^It should be noted, however, that canonical quantization of the system (3) for £ = 1/4 as a system with a 
higher derivatives, should account for its degeneracy (hJ. 
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The energy of vacuum could be determined directly, without resort to eq.(9), as a sum of 
half-frequences interpreted e.g. with the help of zeta-regularization method fLl|. At the same 



time, Green function method JI], §, [11]] used to obtain vacuum energy should rely on the tensor 
T^ v (3) |T0| , but not its first three terms 

= - [d^, d v ip} + - ^g^(d(f) 2 + ^g^m 2 ^ 2 . (11) 

11-component of tensor (3) gives Casimir pressure P = (Tu)sj which, because of translational 
symmetry, coincides with (Tii)«5,z (i-e. the divergence term in (3) does not affect the pressure). 
The role of divergence term is evident from the fact that value £ = 1/4 is the only one at which 
equality 

-^fr = {fu)s > 1 (12) 

holds. Notice that conservation of the total energy |K| and equality (12) both rely on the value 
£=1/4- 

In proving the eq.(12) it is essential that vacuum energy per unit area of the boundary 

~ rl~yo ~ ~ ~ ~ 

E vac = / dxi (Too) 6,i = E c (l, 5) + E wl (x , 5) + E w2 (y , 0) H , (13) 

where dots denote terms vanishing at xo, yo — » as well as at I — ► oo. Only finite part of E vac , 
i.e. Casimir energy 



OC 

E c {l,5) = K D J dt 



(t 2 -/x 2 )V e - 
shtl + tchtJ 



(14) 



depends on the size of the domain and vanishes when I — > oo. Boundary divergences are present 
in r.h.s. of eq.(13) in the form of the energies of the walls, E w \^- For example, 



OC 



E wl (x , 5) = K D [i D ~ l I § (e - 1)^ 1-4 e- 2 — (15) 

The specific representation of energies E w \^ depends on regularization employed, but Casimir 
energy is independent of the latter. Notice that the terms discarded in the r.h.s. of eq.(13) 
have the property of "double vanishing" (with respect to limits of x , y and /) so that energies 
E w \,2 of the walls are determined uniquely within the given regularization scheme. 

By a complete analogy with E vac , one can find its bulk (E vac ) and surface (IT mc = E vac —E vac ) 
components, the first dertermined with the help of density Too (H)- For either of these compo- 
nents the expansion of the form (13) exists giving rize to Ec(l, 5) and Hc(l, 5) correspondingly. 
For example, 

°r t 2 (t 2 — n 2 )~^ 

U C {1, 5) = Kq(1 — D) / dt ^ ?U (16) 

V ; V J J (t + l)(shtl+tchtl) V ; 

and E c = E c — Hc- Since Hc(l, 5) ^ 0, derivative —dE c /dl does not coincide with the Casimir 
pressure (12). Nevertheless, 

*-^-n , (17) 
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and, with the notations A = 5 1 , f = Hc/X, one finds that Ec and Ec are interconnected 
through Legendre transformation: 



Ec(l,X)=E c (lJ(l,X)) + W,X), 



so that 



'dE. 



c 



'dE, 



c 



01 



01 



and 



'dE, 



c 



df 



-A. 



(19) 



/ \ / A 

Now energy Ec determines the Casimir pressure but under specific condition of constancy of 
the quantity / being one half the "Casimir" part of (y?(0) 2 ), see (2). 

3. Asymptotic properties 

For the purpose of comparison with electrodynamics, here we consider a massless case. The 
behaviour of integral (14) in Dirichlet (8 <C I) and Neumann (8 ^> I) regimes is displayed by 
the following expansions (m = 0): 



E c (l,6) = K D 



{2l) D - 1 



-l 



r 2 r 2 + r 3 



1 + 



CP + 2) 
2CP) 



-3 



E C {1, S) = K D t 



-D+l 



A D -{D- l)A D - 2 l, 



D 



3,4, 



D = 2. 



(20) 



(21) 



A 2 + I\n(4r/ E IM - I, 

Here C™ denotes binomial coefficient, A D = 2 1 -°(1 - 2 1 - D )T(D)((D), lnj E = 0.577, and 
eq.(21) contains the leading corrections only. The expansion of pressure (Tu)s } i corresponding 
to (20) has the form (8 < /): 



(Tu) 



5,1 



(D — 1)T(D)((D) 



{21) 



-D . 



(22) 



chi- 1 



CP + 2) 



l- 3 + 



and for Neumann regime (21) 



(Tu) 



s.i 



(D - 1) I 



-D 



2(47r)^r 



D+l 



A D -{D- 2)A D _ 2 I + 



(23) 



The four terms in the square brackets of eq.(22) (term with ^-functions excluding) could be 
obtained through the Taylor series for (I + 5)~ D . Thus, eq.(22) demonstrates the role of 8 as 
penetration depth |]]. Unlike Dirichlet case (22), the correction term in the brackets of eq.(23) 
emerges due to surface energy Uq only. Numerical analysis of the formulas (14) and (16) taken 
at m = shows the dominance of the surface contribution lie over the bulk one (Ec) in the 
region I ~ (1 -=-4)5. The difference in signs between those quantities is responsible for the shift 
of the position of the minimum from I ~ 4.58 (for Ec) to I ~ 1.85 (for Ec)- At that time, 
the qualitative behaviours of energies Ec and Ec as functions of I are alike showing a typical 
van-der-vaals character. 
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At D = 4 formula (22) takes the form 



(Tn) 



Ti' 



5.1 



480Z 4 



4- 



.6 
I 



5 2 ( 
10- + 20 1 



9tt 2 \ 5 3 



185/ I 3 



+ 



(24) 



and should be compared with the corresponding expression in the case of electromagnetic field 
confined to the region between impedance walls (5- skin-depth parameter): 



7T 



2M 4 



165 5 2 640 / 9tt 2 \ S 3 

1 + TJ- 24 F + - 1 + ^'- 4 



740/ I 3 



(25) 



Formula (25) except the last term ~ 5 3 /l 3 in the square brackets, was extracted from [|T], p. 
The corresponding coefficients of expansions (25) and (24) are related to each other in ratios 1, 
1.33, 2.4, 3.46, hence showing the growing effect of spin when going deeper into the boundary. 

4. Interpretation of surface singularities at 5 = and oo 

The surface singularities of the vacuum EMT is a stumbling-block problem for any field theory 
expected to establish the connection between its local properties and observational predictions 
||. It is shown below that the method of dimensional regularization could be applied not only 
to interpret the singular sums like E vac = but gives a reasonable finite answer for the 

V 

energy density of vacuum as well. The method works for even D,5 — +0 or oo and leads to 
an energy expression consistent with the one of ^-regularization method. 

Boundary divergences are represented in r.h.s. of eq.(13) by the energies of the walls (E w i^)- 
Correspondingly, we consider 00-component of a tensor (9) in the limit I —>■ oo, 5 = +0 (alter- 
native case 5 = oo differs from (26) only in sign): 



m 
2n 



R f p- D ' 2 K D/2 (p), 0<p<l, x/l 



0: 



(p')- D/2 K D/2 {p'), <//<!, x'/l^O 



(26a, b) 



(p = 2mx,p' = 2mx' = 2m(l — x)). Exploiting recurrence relations between modified Bessel 
functions K v {p) |TB[, let us transform energy density (26a): 



D 



d 



(^Oo)o.oo "~ > (^0o)o,oo — ~dp^ D ^ ~ ( 27T 



P 1 ' 6 ^e-l(p) 



[l-D)(3-D)-..(l-2e) 



(27) 



Here e = y 



fn( P ) 



n 

2^ j 



0, n — 1, 2 and function 



D 



P 2 %Kn x (p) 



+ . . . + 



P l - £ K e (p) 



;i -D)(3 -£>)... (l-2e) 



l-D (1-D)(3-D) 

possesses following properties: i) for 5fe.D < (3?e e < —n) 

/d(0) = f D (oo) = 0, 
4 ^The singular nature of Dirichlet limit (contrary to Neumann one) is explained in 



(2f 



(29) 
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so that divergence addition in the l.h.s. of eq.(27) does not affect the vacuum energy of the 
wall § 0, 



E wl (0,0)= — / dx p- D ' 2 K D/2 (p) = — — _y^. (30) 
\27i J Jo 8(4vr) — 

(the vacuum energy of Neumann wall E w \(0, oo) = —E wl (0, 0) = —E wl (0, 0)); 

ii) a term f D (p) in the l.h.s. of eq.(27) taken at "physical" values of D — 2, 4, 6, • • •, acts as 
a counterterm eliminating all nonintegrable singularities of density (26a); 

in) this counterterm preserves an exponential decreasing behaviour of the modified density 
(27) at x > %/mc. The latter is not still uniquely defined: one can add to /d(p) any regular 
function having the property (29) (see || 

Thus, at 5 = +0 the total energy of vacuum includes Casimir part |14 



m D l 7^(e 2 -l)^ 
(47r) £ f i r 



E a(h 0) = "77-^7777177 / « = 1 (31) 



and doubled topological (according to 0) energy (30) caused by the presence of the edge of 
the manifold. At 5 = oo the total energy of vacuum reduces to Ec(l, oo) because the energies 
of Neumann (x = 0) and Dirichlet (x = I) walls cancell out. Energy Ec(l, oo) = Ec(l, oo). 
The integral representation for it can be obtained from (31) by substitution "+" for "— " in the 
denominator of the integrand and by multiplying the whole expression by (—1). 

Dimensional regularization method does not result in satisfactory expressions for the vacuum 
energy when < 5 < oo. As it may be seen from eq.(15) (taken at xq — 0), for even D the 
corrections of order p, p 2 , • • ■ to the energy (30) cannot be interpreted in the way like (30). On 
the other hand, for odd D > 3 the energy (30) is infinite. Now, the depending on p " corrections" 
find finite interpretation, but being summed up, they lead to logarithmic divergence in the 
Neumann limit p —>■ oo. 

The analytical structure of the surface singularities of the vacuum EMT (26) with respect to 
dimension D is just that property which makes dimensional regularization applicable to define 
finite density (27). In view of the absence of translational symmetry it seems quite natural 
for density like (27) to exist. At that time finite renormalization of the type suggested in 
H [14| eliminates the energy of the walls from E vac = E c (l, 0) + 2E W . This procedure needs in 
explaination how to interpret the expressions like (27). It should be mentioned in addition that 
the energy density (27) (but not (26)) vanishes when m — 0. This disappearance was assumed 



to be characteristic of conformally invariant models pAj |15 



5. Conclusion 

Below some comments on literature related to the present topic are given. The massless 
case for the scalar model at hand (dimension D = 2, 3) was considered in works |], IB] and 



formulas consistent with (9) ( D — 3, m = 0, Hl6| ) and (14) ( D — 2,m — 0, 0) were derived. 
The work |IBJ exploits an idea |T7| of replacement of the boundary with a singular non-locally 
regularized potential, that occured to be equivalent to taking potential energy (2) into account. 
Formula (2.18) from [14[] corresponds to our formula (31). Doubled energy (30) appears in 



TA\ as well, but it was not associated there with the vacuum energy of half-space. Our formula 
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(9) at 5 = +0 (or, equivalently, (31) and (12)) gives the expression 



(2ii>o,i 



m D (l - D) 



(4tT 



e(e - 1) 

e 2mlS, _ \ 



D-3 
2 



(32) 



m 



D 



-2ml 



(A-Kml) 



D-l 
2 



1 | (D-l)(D + 5) 
16m/ 



+ 



which does not coincide with formula (2.13) from The latter, being taken at zero temper 



ature, includes physically unacceptable dependence of pressure on (arbitrary) renormalization 
parameter. Such dependence, as it was noted above, is characteristic of curved boundary under 
the condition of coincidence between the parameters determining its curvature and size. 

The author would like to acknowledge useful discussions with A.I. Nikishov and V.I. Ritus 
and financial support from RFFR (grants 95-02-04219-a and 96-15-96463). 
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